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We experimentally demonstrate a simple method to measure the biphoton joint spectrum by mapping the 
spectral information onto the temporal domain using a dispersive medium. Various top-hat spectral filters are 
used to limit the spectral (and hence, temporal) extent of the broadband downconversion photons measured. 
The sharp edges of the spectral filters are utilized as spectral markers for dispersion characterization of the 
dispersive medium. This method allows dispersion characterization and joint spectral measurement to be 
completed simultaneously. The joint spectrum (which extends beyond 100 nm, centered about 1.5 gm) of the 
type-II downconverted photon pairs generated from a poled optical fiber is obtained with this method. 
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1. Introduction 

Understanding the spectral correlations present in a cor¬ 
related photon pair source is extremely important in 
determining whether that source is suitable for various 
quantum interference and quantum information process¬ 
ing schemes. For instance, a photon pair whose signal 
and idler photons are spectrally-uncorrelated may work 
well as a heralded source M , but be unsuitable for ap¬ 
plications where spectral entanglement is required, such 
as quantum-enhanced optical coherence tomography [|J, 
or high-dimensional quantum key distribution [3, [6[ . 

However, traditional measurements of the joint spec¬ 
trum of a biphoton require the use of expensive and 
lossy spectrometers, involving long integration times to 
produce significant coincidence counts. The integration 
times are increased when the resolution required is finer. 
Recent progress in so-called ‘stimulated’ emission to¬ 
mography [7| has shown that the joint spectrum can 
be predicted accurately by performing classical nonlin¬ 
ear frequency-mixing. In practice [§j, though, this form 
of tomography requires complete knowledge of and ac¬ 
cess to the nonlinear generation process (which will not 
be the case if the source is embedded or packaged in a 
commercial system). 

Another way to approach the problem is to map the 
spectral information of the JSI onto the temporal do¬ 
main, using a dispersive medium, so that the spectral 
information is now encoded in the arrival times of the 


signal and idler photons at their respective detectors. 
When the dispersive medium is a long spool of fiber, 
this is called the fiber spectrometer method [§). How¬ 
ever, the fiber spectrometer requires accurate dispersion 
calibration that is specific to the actual fiber spool in 
use. In HO, for example, the fiber dispersion is cali¬ 
brated beforehand with a short-pulse tunable laser. This 
must be done before every experiment, as the temper¬ 
ature changes in Corning SMF-28 fiber, say, can result 
in a change of group delay by 40 ps/km/K [Hj, or 4 ns 
over 20 km of fiber for a 5 K temperature variation. 

In this work, we demonstrate experimentally that such 
a priori calibration is unnecessary, especially if we wish 
to measure a biphoton joint spectrum that is broadband. 
By using multiple frequency-conjugate filters that have 
sharp cut-on and cut-off wavelengths, the broadband 
spectrum we obtain a new way to elucidate the disper¬ 
sion of the fiber spool. This is because the sharp edges 
of the filter can now act as spectral markers, allowing 
us to simultaneously gauge the fiber spool’s dispersion 
while obtaining the joint spectrum of the biphoton’s JSI. 
Additionally, the use of the filters limits the temporal ex¬ 
tent of the downconverted biphotons, so that a relatively 
high-repetition rate pump (Ti:Sapphire laser, 81.6 MHz, 
12.25 ns period) can be used without the downconverted 
biphotons spanning more than one period, which could 
cause complications experimentally. 

The breakdown of this paper is as flows. In Section 
2, the idea of the fiber spectrometer will be explained 
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mathematically. Using downconverted photon pairs gen¬ 
erated from a poled optical fiber jl2l - fl5| and various top- 
hat sharp-edge filters, we will demonstrate our method 
to measure the dispersion of the fiber spool in Section 
3. Finally, in Section 4, the fiber spectrometer method 
will be used to measure the downconversion spectrum 
of the poled optical fiber. Conveniently, the downcon¬ 
verted photon pairs generated have wavelengths sym¬ 
metric about ~1550 nm. 

2. The fiber spectrometer 

Consider a biphoton generated from SPDC. Both pho¬ 
tons are created at the same time and will also exit the 
nonlinear medium simultaneously; in this case, ‘simul¬ 
taneous’ shall mean within the detector response time 
(120 ps full width at half maximum [FWHM]). Let us 
describe this biphoton as a function of wavelength: 

Nph — / d\ s d\i fi(As,Aj), (1) 


where N p h is the average number of photon pairs in 
the wavepacket, and h(X s ,Xi) is the spectral brightness 
(photon pairs per unit wavelength of signal and idler, 
[pairs/nm 2 ]). 

Now, let this biphoton propagate through a fiber spool 
of length L whose dispersion is given by the D-parameter 


D{ A) = 


2-7TC d 2 jd 


( 2 ) 


where c is the speed of light in vacuum and /3 is the 
fiber propagation constant. The fiber is single-mode at 
the wavelengths of the downconverted photon pairs. We 
now exploit the fact that single photons travel at the 
group velocity [l7] to develop the fiber spectrometer idea 
further. 

The spectral information of the biphoton will then be 
mapped onto the temporal domain: 


N ph = j dt s dti n(t s (X s ),ti(Xi)), (3) 

where we have explicitly stated that the arrival time t s 
( ti ) of the signal (idler) is a function of its wavelength 
A s (Aj). It is a simple matter, then, to relate n (Eq. [3]) 
to n (Eq. |TJ). Noting that the arrival time is defined as 

t s ( A s ) = —7— - (where v q is the group velocity inside 
Vg(X s ) 

the fiber spool) and that 

dt a = d ( v ° {Xs) ^ dX s = LD(X S ) dX s , W 

CLAs 

we have the equality: 

n(X s ,Xi) = L 2 D(X s )D(Xi ) n(t s ,U). (5) 

In our case, this spectral-to-temporal mapping is one-to- 
one because the downconverted photons generated from 


the poled fiber have wavelengths that lie only in the 
anomalously dispersive region of the fiber spool used 
(SMF-28). 

Note that while D{X) is effectively the group-velocity 
dispersion (GVD) of the fiber spool, it encapsulates 
higher-order dispersion terms of the spool as well. 

Experimentally, an actively-modelocked Ti:Sapphire 
laser with an 81.6 MHz pulsetrain and 400-ps FWHM 
pulses pumps a periodically-poled fiber to generate an 
ensemble of biplrotons. Fig. |T| shows an outline of the 
experimental setup. The signal photons (shorter wave¬ 
length) are separated from the idler photons (longer 
wavelength) using coarse wavelength-division demulti¬ 
plexers (DEMUX); the bandwidths of such coarse filters 
range from 16 nm to 40 nm (see Fig. E|). Free-running 
SPDs, the id Quantique id220 devices, are used to mea¬ 
sure the photons. The arrival times of the signal and 
idler photons at the SPDs are measured with respect 
to the sync signal of the Ti:Sapph with the commercial 
time-interval analyzer (PicoQuant HydraHarp 400). A 
two-dimensional histogram is constructed from the de¬ 
tector clicks and Ti:Sapph sync signal. A schematic for 
the process is shown in Fig. [3] 

Equation fallows us to relate the temporal resolution 
to the differential wavelength resolution dA m i n . The min¬ 
imal resolution dt min is dictated by the jitter of the sync 
pulse (cr sy nc), the pump pulsewidth (cr pu ise), as well as 
the jitter of the detectors (troet)- We make the (good) 
approximation that all three are normally-distributed 

exp 2 ^)), so that the jitters will add in quadra¬ 
ture to give the minimal resolution dt m i n : 


dtmin — 


pulse 


+ °Def 


( 6 ) 


where <7 syn c < 30 ps, cr pu i se = 170 ps (as FWHM ss 
2.355cx), and UDet = 51 ps. We see that the expression 
is dominated by er pu ise- The value of dtmin (180 ps) 
corresponds to a wavelength resolution dA m i n of 0.13 nm 
at 1550 nm for a 20-krn spool of Corning SMF-28 fiber. 

Finite-bandwidth filters are required as the range of 
the arrival times of the signal photon must be no greater 
than 24.5 ns (or one period of the laser sync signal, 40.8 
MHz). Multiple filter sets (Fig. [2]) allow us to span 
a large portion of the downconversion spectrum (more 
than 100 nm). 

An average photon pair N p h of no more than 0.01 
pairs/pulse is generated over each filter set’s bandwidth. 

Using Eq. 0 we obtain the ‘raw’ spectral brightness 
nRAw(X s ,Xi ); we say it is raw because we still need to 
account for the spectral dependence of the insertion loss 
of the system (including detector efficiency and the sys¬ 
tem insertion loss), and normalize by the average pump 
power and detector deadtimes in order to obtain the 
true spectrum of the biphoton. Finally, after obtaining 
n( A s , A,), we can relate it to the biphoton wavefunction 
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Fig. 1. A simple schematic of the fiber spectrometer setup. 
The spectral information of the signal and idler photons are 
mapped onto the time domain after traveling through the 
long fiber spool. A WDM (labeled DEMUX) is used to route 
the shorter-wavelength signal to one detector, and idler to 
the other detector, while also limiting the spectral extent 
of each photon. This is important, as it limits the temporal 
extent of the signal photon, which cannot exceed 24.5 ns (the 
sync period of the Ti:Sapph oscillator). 

4>(X s ,Xi), where: 

\il>) = jdx s dXi 4>(X s ,Xi), (7) 

in the following way: 

h(X s ,Xi) = ji\<f>(X s , Ai)| 2 . (8) 

The factor jl is the total brightness, defined as the av¬ 
erage number of photon pairs generated (over the entire 
downconversion spectrum) per mW of pump power. 

3. Fiber Dispersion Characterization 

The dispersive medium used here for the spectrometer is 
a 20-km reel of SMF-28 fiber. While the GVD of SMF-28 
is well documented, there may be some manufacturing 
variability from one reel to another. In this subsection, 
we detail how we measure the dispersion of a fiber spool 
used for the experiment. But first, we must assume a 
form for the D parameter 0: 

°W = f (* - I) (9) 

where Ao is the zero-dispersion wavelength, and So is the 
dispersion slope at Ao- 

We wish to know So and Ao for our particular fiber 
spool, as well as the suitability of Eq. [9] over the wave¬ 
length range of interest (1480-1620 nm). To that end, 
we can exploit the filtered downconverted light from our 
biphoton source to characterize the dispersion of the 
fiber spool. We observe that in Figs. ED: and[3Ji, the one¬ 
dimensional histograms of the signal and idler arrival 
times (respectively) have sharp rising and falling edges 
(> 15 dB/nm). These sharp features correspond to the 
cut-on (Ai) and cut-off (A 2 ) wavelengths of our filters, 
and using multiple filter sets (1490-1610, C-/L-Band fil¬ 
ters, 1530-1570) to filter the downconverted light, we can 


perform a chi-square fit, using the formula below (Eqn. 

CUD: 


At = L [ dX D (A), (10) 

J Ai 

with the cut-on/off wavelengths as the independent vari¬ 
able, and the range At of arrival times as the depen¬ 
dent variable. The best fit yields the following results: 
S 0 = 0.0885 ±0.0006 ps/(nm 2 -km), A 0 = 1320 ± 1.7 nm, 
and L = 20.56 ±0.13 km. 

With the fiber spool dispersion characterized, we 
are ready to use it to invert the temporal histogram 
data n(t s ,ti ) to obtain the biphoton spectral density 
n(X s ,Xi). 

4. Experiment and Results 

The fiber spectrometer method is used in the measure¬ 
ment of the joint spectrum of a broadband correlated 
photon pair source, namely the poled optical fiber E2- 
Il5p . The poled fiber has a non-zero second-order nonlin¬ 
earity, and is quasi-phase-matched [0 for the downcon¬ 
version of biphotons at ~ 1.55 pm. The type-II down- 
conversion process [10 is exploited; however, neither 
the polarization-entangled nature nnnim nor the 
polarization degree-of-freedom in general are exploited. 

The source’s phase-matching peak occurs at the pump 
wavelength X p o = 774.5 nm, with a FWHM pump accep¬ 
tance bandwidth of 0.35 nm. The pump pulsetrain used 
is narrowband (< 0.02 nm) compared to this bandwidth. 
Due to the broadband nature of the downconverted pho¬ 
tons, and the narrowband pump used, the generated 
biphotons are highly spectrally-correlated (Fig. [2}i). 
Additionally, the biphoton spectrum varies significantly 
as a function of pump wavelength. 

The experimental procedure is as follows. The pump 
is set to a particular wavelength A p , and generates down- 
converted biphotons when launched into the source. The 
biphotons are then coupled into the dispersive fiber spool 
(Fig. HD, and then filtered with frequency-conjugate fil¬ 
ters (Fig. HD at the output. A 2-D temporal histogram 
(n(t s ,ti)) tallying up the coincidences as a function of 
signal and idler photon arrival time is collected (Fig. 
Eb)- Integration times of 500-1000 seconds are used to 
generate each histogram. The temporal extent for the 
2-D histogram is measured along the t s (A t s , Fig. EDO 
and t,j axes (A t ,, Fig. EJi), as these values are used to 
obtain the fiber spool’s dispersion characteristics (Eqn. 
HUD. The procedure is repeated with another pair of fil¬ 
ters or pump wavelength until a 120 nm range at the sig¬ 
nal/idler wavelength and a 1 nm range about the central 
pump wavelength A p o is swept (see Fig. HDi). The pump 
wavelength is varied at a step size of approximately 0.1 
nm. 

At certain pump wavelengths, not every filter pair 
will be completely frequency-conjugate. For example, 
consider the C-Band (L-Band) filter, that has cut-on 
and cut-off wavelength from 1498 to 1540 nm (1567.6 
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y = 1530±8,1570±8 nm 
8 =1528-1563 nm 



Fig. 2. a) The various top-hat filters used are overlayed on top of the calculated joint spectra for the downconverted photon 
pairs, b) The (relative) transmission spectra of the filters used are plotted, a) X s = 1482 — 1498 nm, A i = 1602 — 1618 nm. /3) 
As = 1498 — 1540 nm, A i = 1567 — 1612 nm. 7 ) A s = 1522 — 1538 nm, Ai = 1562 — 1578 nm. 5) The degenerate filter; the filter 
is from 1528 nm to 1563 nm. 


to 1612.6 nm). Tabled] gives the signal filter’s actual 
frequency-conjugate bandwidths (from Ai to A 2 ) and 
temporal extent A t s (Fig. [3J:) at various pump wave¬ 
lengths. 


A p 

Ai 

A 2 

Time-Width 

(nm) 

(nm) 

(nm) 

A r 3 (ns) 

774.7 

1498 

1531.7 

9.77 

774.9 

1498 

1532.5 

10.05 

775.1 

1498 

1533.1 

10.20 


Table 1. Data obtained using the C-Band filter as a function 
of pump wavelength. The cut-on and cut-off wavelengths 
refer to the portion of the filter where photons that have 
conjugates residing in the transmission window of the L-Band 
filter are generated. The overall bandwidth for the C-Band 
(L-Band filter) is 1498-1540 nm (1567.6-1612.6 nm). The 
error for A t 3 is taken to be di m i n (see Eqn. or ±180 ps. 

Data similar to Table[l]is collected from all the experi¬ 
mental runs for various pump wavelengths and filter sets, 
and fed into Eqns. l9land [T0l from which the dispersion 
properties (D(A)) of the fiber spool are then determined 
(as described in Section 3). 

Once the D— parameter of the dispersive spool is 
known, the 2-D temporal histogram n(t s ,ti ) (see Fig. 
[3}o) can then be related via Eqn. [5] to n(X s ,Xi). The 
data for the various filter sets and pump wavelengths is 
then stitched together to give a tuning curve (Fig. 0}i). 
The x— axis is the pump detuning A (A = X p — A p o), 
while the y— axis is the signal/idler wavelength. A cross- 
section of this tuning curve (at A = 0.25 nm) is shown 
in Fig. HJo. 

In summary, we have demonstrated a method for mea¬ 
suring the broadband downconversion spectrum of a cor¬ 
related photon pair source using a fiber spectrometer. 


The dispersion of the fiber spool can also be simulta¬ 
neously measured using the same source of broadband 
correlated photon pairs. 
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